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We investigate a classical phase-space approach of matter-wave propagation based on the Trun- 
cated Wigner Equation (TWE). We show that such description is suitable for ideal matter waves 
in quadratic time-dependent confinement as well as for harmonically trapped Bose Einstein con- 
densates in the Thomas-Fermi regime. In arbitrary interacting regimes, the TWE combined with 
the moment method yields the low-energy spectrum of a condensate as predicted by independent 
variational methods. TWE also gives the right breathing mode frequency for long-ranged interac- 
tions decaying as 1/r 2 in 3D and for a contact potential in 2D. Quantum signatures, beyond the 
TWE, may only be found in the condensate dynamics beyond the regimes of classical phase-space 
propagation identified here. 
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Introduced in the early days of quantum mechanics, 
Wigner distributions have been applied successfully since 
then in many areas of physics. In optics it is used to de- 
scribe partially coherent light beams [l|, and in atomic 
physics to characterize atom interferometers [2j or Bose- 
Einstein condensates (BECs) 3]. The Wigner distribu- 
tion obeys a propagation equation involving a series of 
differential operators weighted by increasing powers of 
h [J] which can be identified as the quantum terms of the 
Wigner equation. The classical limit of this equation, 
obtained by taking the limit h — >• 0, is referred to as the 
Truncated Wigner Equation (TWE). This equation has 
been used with stochastic classical fields to study quan- 
tum gases @ . We propose here to use it instead with the 
modes of a coherent atomic beam. When combined with 
the moment method, this reveals an accurate yet simple 
way to capture the low-energy dynamics of a BEC. It 
is for instance remarkable that the TWE describes cor- 
rectly the time-of-flight expansion of condensates in the 
limit of vanishing interactions as well as in the Thomas- 
Fermi limit. Our treatment, which studies the agree- 
ment between the TWE and the full Wigner equation, 
also identifies the circumstances under which quantum 
signatures -i.e. effects due to the quantum terms of the 
Wigner equation - may appear in the dynamics of BECs. 

After a quick reminder on the general Wigner equa- 
tion, we study its applications to the propagation of ideal 
matter waves in time-dependent quadratic potentials. 
We retrieve the ABCD propagation formalism [6] for 
coherent guided atom optics [7|, suitable to investi- 
gate atomic beams propagating in time domain [8[ 
through a variety of atom-optical elements. It can 
also be used to study the transverse stability [9[ of an 
atomic resonator [lOj , to define the quality factor of an 
atomic beam [ll|, or to investigate the generation of 
atom-optical caustics [ljj caused by a sudden potential 
change. We also analyse the predictions of the TWE 
for the expansion and for the low-energy spectrum of 
interacting atomic waves evolving in harmonic traps. 



The consistency of the TWE with a universal prediction 
for long-ranged 1/r 2 interactions is verified. 

We recall the form of the Wigner Equation in the pres- 
ence of a general two-body potential V(r). The non- 
relativistic atomic field operator ip obeys the general 
propagation equation 

iK-fi(r,t) = -— A^(r,t) + J7 (r,t)^(r,t) 

where Uo(r,t) is the time-dependent external potential 
experienced by the atoms. We assume the sample to be 
well-described by the mean-field approximation and in a 
number state associated with a macroscopic single mode 
<j>. The associated Wigner distribution is 



W(r,p,t) 



dr' 
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The mode <f> satisfies a Schrodingcr equation with an 
effective Hamiltonian H(r,p,t) = p 2 /2m + U(r,t). 
where U(r,t) = U (r,t) + C/ mf (r,t) and U ml (r,t) = 
J dr'\cf>(r' ,t)\ 2 V(r — r'). For a contact interaction po- 
tential V(r) = gS(r), this Schrodinger equation re- 
duces to the usual time-dependent Gross-Pitaevskii equa- 
tion(GPE). The Wigner distribution associated with a 
wave- function governed by an Hamiltonian H(r, p, t) sat- 
isfies the transport equation [12|, [l3| dW/dt = —i£ [W], 
where C is the Liouvillian operator 



C = H(r,p,t) 
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We have identified the vectors r — (x, y, z) and p = 
(Px,P y ,Pz) with (ri,r 2 ,r 3 ) and (pi,p2,p 3 )- The differ- 



entiation operators a and d act on the Wigner dis- 
tribution W(t, p, t) and on the Hamiltonian H(r,p,i) 
respectively. The Wigner equation can be recast in an 
explicit differential series ordered by increasing powers 
of the constant h 
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the following ansatz as a possible solution to the TWE 

F(r,p,t) = l^( T ^(r-t)-- T S(p-™^), 
V m 

- m T C(r - ~i) + T A(p - m"?) , to) (7) 

This ansatz indeed satisfies the Wigner equation if the 
following relation is fulfilled for any phase-space point 

Tir • W° + T 2 — • W r ° + mT 3 r • W° 
m v 

+ T 4 p • W" + f 5 ■ W? + mf 6 ■ W° = . 



To alleviate notations, from now on we note the 
partial differentiation with respect to time, space 
and momentum coordinates thanks to an index (i.e. 
W 



J p 2 



a" 



dp 1 dp 2 dp 

classical limit h 



is the derivative of the Wigner distribution 
W). 



The TWE is obtained by taking the 

0, i.e. by discarding the last term of 
the right-hand-side (r.h.s.) in Eq.(j4|). 

Let us first treat the propagation of ideal atomic waves 
in time-dependent quadratic potentials. In this regime, 
the full Wigner equation exhibits no quantum terms and 
thus reduces to the TWE. We show that the Wigner 
equation, or here equivalently the TWE, can be used 
to recover exactly the ABCD propagation formalism. 
This phase-space method had been already verified by 
direct integration for ID Laguerre-Gaussian modes 
and analogously for partially coherent matter waves 



14], 
15] 



The Wigner equation offers a compact proof of this re- 
sult, avoiding the computation of propagation integrals. 
Consider a general quadratic Hamiltonian (without 
atomic interactions) expressed with the conventions of [6| 

T p/3p T m T . -i 

H = — rap - — r 7 r - mj • r + / • p , 5) 

2m 2 

a, ft and 7 are 3x3 time-dependent real matrices sat- 
isfying T a — —a, T (3 — (3, T/ -f — 7, and f ,g are time- 
dependent vectors. T stands for the matrix or for the vec- 
tor transposition. One sees readily that the full Wigner 
equation is here simply given by its classical part (TWE) 
Vt = —Hp( r ,p) ■ 2/r + H r (r,p) ■ y p . We introduce the 
phase-space vector R(i) = (r, p/m). We expect the fol- 
lowing phase space map @ 



R(t)= M(t,to)K(t ) + R {t,t ) 



(6) 



Ro(£)£o) = ( £ 5 4> ) is a source term to be determined 
later. M(t, to) is a 6 x 6 matrix, referred to as the ABCD 
matrix associated with the evolution between the initial 
time t and final time t: M(t,t ) = [A\B,C\D] . We 
used the matrix notation M = [mii|mi2,m2i|m22] - m ij 
being the element of the i^ row and j^ column- and 
A, B, C, D are 3x3 matrices depending on the pair of 
instants (t, to)- The ABCD matrix satisfies the symplec- 
tic relation M~ x = [ T D\ - T B, T C\ T A]. The phase-space 
map ^ and the symplectic relation lead us to consider 



with 

Ti 
T 3 

J 5 
-± 

T 6 



= 'D 

= - T C 



T D T a- 
T C T o 



l B~f, T 2 = - 1 B + 1 D 1 P 
- T A 1 , T 4 = T A- T Cf3 - 



[ Ba 



T A T o 



l B, 



<-> 



,-> 



l Bcf,- T Bt-{ T Dt- T D-t + ^D'j) 

= T C~t + T ct - T c7 - ( T i t + T A~i T A~f) 



The upper dots stand for time derivatives and lower dots 
denote scalar products. The terms W®, W® correspond 
to the gradients of the Wigner distributions towards the 
position and momentum vectors respectively, evaluated 
at the initial time to. It is now convenient to introduce 
the matrix N and vector S defined from the Hamilto- 
nian © as N = [a|/3,7|a] and S = (/ \~$). Eq.© im- 
plies that Jj. 2,3,4 = O3 and T 5,6 — 1 i- e - the parameters 
A,B,C,D, £ , and <j> satisfy the differential equations 



AM 



NM, -(M- l R ) 



M~ l S, 



(8) 



with the initial conditions £, (t 0l t ) = 4>(t ,t ) = 0, 
A(t ,t ) = D(t ,t ) = h and B{t 0l t ) = C(t ,t ) = 3 . 
I3 and O3 denote respectively the identity and the zero 
3x3 matrix. The system obtained in [6] for the pa- 
rameters A, B,C,D, £ , (j> is equivalent to (jHJl. With the 
adequate initial values for A,B,C,D, £ , cf> , the ansatz F 
coincides initially with the Wigner distribution and sat- 
isfies the same differential equation: the Wigner function 
thus evolves according to Eqs. (J7J and ©. 

The TWE also gives the right predictions for the ex- 
pansion of a dilute BEC trapped in the Thomas-Fermi 
limit, in particular the correct expansion law |16| ex- 
pected in a time-of-flight experiment and the correspond- 
ing low- lying mode frequencies |17| . Here, we summarize 
only the main points of the argument. The BEC, sub- 
ject to contact interactions, evolves in a time-dependent 
quadratic external potential U(r,t) = J2i=i TOW i 2 (*) r ?/2- 
Our approach combines the technique of scaling fac- 
tors already exploited by several authors [16|, Oil Qjl 
with the TWE. The transposition of scaling laws to 
Wigner functions consists in searching for an ansatz of 
the form W(r,p,t) ~ W(r' ,p',0) with a gauge trans- 
form 20] (rj,p{) = {ri/bi{t),bi{t)pi - mbi{t)n). This 



ansatz depends on time only through the scaling param- 
eters hi. Its adequacy with the exact Wigner distribution 
in the Thomas-Fermi limit is not obvious and has not 
been considered so far to our knowledge. For this pur- 
pose, one writes the Wigner distribution as W(r,p,t) = 
W(r',p',0) + SW, and shows that the time-dependent 
correction SW vanishes in this regime. Inserting this 
expression in the TWE, one sees that this correction is 
driven by a term which vanishes if the scaling parame- 
ters satisfy b\ = u>f(0)/(6, Tlj=i fy) — &jW 2 (t). Using a 

procedure suggested in Ref. [Iff, one can show that SW 
remains negligible in the Thomas-Fermi limit. 

A powerful theoretical technique to track the phase- 
space motion is provided by the use of moments, and 
has been extensively used in combination with non- 
linear Schrodinger-like equations, in optics and electro- 
magnetism [2l|, [22| . It turns out to be useful also to de- 
rive low-lying collective modes of a harmonically trapped 
Bose gas degenerate or not. In the following, we detail 
the transposition of this technique within the Wigner's 
formalism. Moments are defined with the Wigner distri- 
bution by (o(r,p))(i) = J d 3 rd 3 p o(r, p) W(r,p,t) with 
a(r,p) a generic monomial in the position and momen- 
tum coordinates. By convention we define the order of 
the moment as the total power of the monomial. The 
orders in position and momentum are the partial powers 
associated respectively with the product of the position 
and the momentum coordinates inside the moment. 

It is insightful to investigate whether the predictions of 
the TWE regarding the dynamics of moments differ from 
those of the full Wigner equation. From the full Wigner 
equation Q, we deduce that the first-order time deriva- 
tive of moments with an order in momentum strictly less 
than three can be expressed exactly by considering only 
the classical terms of this equation, that is the TWE. In- 
deed, considering the moment (r^r^r^p^p^p^ 3 ), one 
sees that its time derivative involves only the terms of the 
quantum series ^ satisfying nu < bk for k = 1,2,3. In 
particular, the first-order time derivative of any second- 
order moment is given exactly by the TWE. 

For ideal matter waves propagating in quadratic po- 
tentials, the two terms in the r.h.s. of the TWE, namely 
—pi/m W ri and U ri W Pi , couple the motion of a moment 
only to other moments of equal or lesser order. With the 
remark above on the time derivative of moments, one 
obtains that the TWE describes exactly the dynamics 
of second-order moments. Besides, any moment follows 
also a closed set of equations of motion. In contrast, if 
the effective potential U contains polynomials beyond the 
second order - coming either from the external potential, 
either from mean- field interactions -, these properties are 
no longer verified. Indeed, the second classical term in 
the r.h.s. - U ri W Pi - couples the motion of any moment 
with a non-zero momentum order to higher-order mo- 
ments. Second-order moments are then coupled to a hier- 
archy of moments of arbitrary high order, some of which 
influenced by the quantum terms of the Wigner equa- 



tion. In this case, one must use the full Wigner equation 
to capture the moment dynamics. 

By applying the TWE to specific moments which are 
not directly influenced by the quantum terms, it is still 
possible to infer exactly the low-energy mode frequencies 
of an interacting sample in the presence of either contact, 
dipolar or 1/r 2 interactions. We focus from now on the 
motion of second-order moments. Their first-order time 
derivative can be evaluated exactly without considering 
the quantum terms in the Wigner equation. In this sense, 
the equations of motion [23| 



dt 



(npi 



d(riPi) 



1 



dt 



(PJ 



■mwAr;) 



■<riU%*) 

(9) 
are classical and can be regarded as TWE predictions. 

Let us prove that the TWE respects the universal 
ground frequency invariance 24] expected when one adds 
an interaction potential V satisfying the scaling law 
V(Ar) = A~ 2 U(r). In a 3D system this property is only 
fulfilled by the long-ranged 1/r 2 potential, but in a 2D 
system (such as a condensate with a frozen external de- 
gree of freedom) it is also valid for contact potentials such 
as gS(r). Similarly to the derivation of the virial theorem 
in mechanics, one obtains by differentiation of Eqs.(|9]) 



m d 2 (r 2 ) (p 2 ) 2 2 

2 dt 2 m ov ' 



U, 



mf\ 



(10) 



The scaling law of the interaction potential implies that 
(r • C/ r mf ) = -2(f/ mf ). One obtains the closed equa- 
tion d 2 (r 2 )/dt 2 + 4w 2 (r 2 ) = 4E/m with the total energy 
E = (p 2 )/2m + muj 2 (r 2 )/2 + (U ml ). The lowest mode 
frequency 2wo is thus immune to the addition of an in- 
teracting potential with the scaling property above. 

We now show that the classical Eqs.© are sufficient to 
predict the low-energy modes of an interacting conden- 
sate. We choose an interaction potential in the general 
form [25[ V(r) = g c 5(r)+g d (l — 3 cos 2 r )/r 2 . It contains 
a contact term and a dipolar contribution associated with 
a permanent moment, dipoles are aligned along the vec- 
tor e3 and cos 8 r = (r-e^/r. As such the TWE cannot be 
solved analytically, but we can investigate its predictions 
thanks to a Gaussian ansatz 



W{v,p,t) 



(nh) 3 11 

V ' »=1,2,. 



Hli) 



-{pi-hai-2hPi Tif 



(11) 



rot is the average position, Xi the width of the wave- 
packet, and the parameters on,fii define a quadratic 
wave-front. Using the relation d{ri)/dt = (pi)/m and 
the first Eq. ©, one finds a, = (m/h)[f Qi — (Xi/Xi)roi\ 
and j3i = mXi/(2HXi). Using these expressions in Eq. fTTj) 
yields the scaling gauge transform introduced previously. 
Eqs.© give for each coordinate an equation of motion 
similar to Eq. (fP0]) . The second-order moments in posi- 
tion and in momentum can be written with the ansatz pa- 
rameters as (r 2 ) = r 2 i + X 2 /2 and (p 2 ) = m 2 (f^ + A 2 /2) + 
h 2 /(2X 2 ). Reporting these expressions in the equations 



of motion, and using roi = — wfroj, one finds that the 
Gaussian widths A satisfy 



mXi = —mui; Xi 



tA? 



A, 



(n EC) 



(12) 



The contact potential contributes to the r.h.s. with 
(nil?.) = ~g c N/(4V2n 3 / 2 X 1 X 2 X 3 ), yielding a gradient 
in Eq. (fl"2|) . By setting ^ = and Aj = i?f&i with 
adequate constants Up, Eq. (fT2"|) coincide with the pre- 
vious differential system for the parameters 6, up to a 
term vanishing in the Thomas-Fermi limit. The Gaus- 
sian ansatz, successful in both the dilute and strongly 
interacting limits, is thus a valid interpolation between 
these regimes. This also shows the resilience of surface 
modes towards the sample profile. The dipolar poten- 
tial also contributes through a gradient, thereby turning 
Eq. (|12p into the equation of motion of a fictive point-like 

particle of mass m, of position A , and experiencing an 
effective potential V e 



ftf)=E( 



h 2 m 2 2 \ N 

V2^A| + Y Wj V + 7t3/ 2 A 1 A 2 A 3 



7! + W ^/ dr 



1 - 3 cos 2 i 



.J] c -r?/(2A?) 



(13) 



This result matches perfectly the predictions of the vari- 
ational methods 19, 26], in which the same differential 

system m Xi = — dV e (X)/dXi was obtained. One 
retrieves the low-energy monopolar and quadrupolar fre- 
quencies for a cylindrical condensate obtained by other 
theoretical methods I19L l27j and confronted successfully 
to the experiments |17| . The method of moments com- 
bined with the TWE is thus entirely equivalent to the 



variational method on this class of Gaussian ansatz, but 
it nicely avoids the usual algebraic operations required 
in the latter to uncouple the Euler-Lagrange equations. 



In conclusion, the classical phase-space evolution given 
by the TWE is sufficient to explain many important fea- 
tures in the dynamics of zero-temperature condensates. 
Using an adequate transformation of the Wigner distri- 
bution, we have recovered the dynamics of an expanding 
BEC in the Thomas-Fermi regime. With the method 
of moments, we have obtained the correct low-energy 
spectrum of a condensate in the presence of contact and 
dipolar interactions of arbitrary strength. Besides, the 
universal oscillation frequency in the presence of a 1/r 2 
potential is also predicted by the TWE: this suggests 
that the hidden symmetry of the full equation (2J| is pre- 
served despite the truncation of the Wigner equation. 
This phase-space approach has allowed us to revisit the 
ABCD law for the propagation of dilute and partially co- 
herent matter- waves, which can be extended to account 
for mean- field interactions [9j. Last, we have shown that 
the quantum terms of the Wigner equation do not affect 
the low-energy modes of trapped interacting condensates. 
Since quantum signatures appear in regimes where the 
TWE fails, the presented results suggest that neither the 
Thomas-Fermi nor the diluted limit are appropriate to 
evidence such signatures: these should be rather tracked 
with intermediate mean-field interactions or in high-order 
modes. Recently, it has been brought to our knowledge 
that the variational and moment methods have been also 
investigated in [281 ] . 
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